The Sunyaev-Zeldovich Signal of the maxBCG SDSS Galaxy Clusters in 

WMAP 



Patrick Draper^, Scott Dodelson^'^'^, Jiangang Hao^, Eduardo Rozo'^''* 
^Enrico Fermi Institute, The University of Chicago, Chicago, IL 60637 
^ Center for Particle Astrophysics, Fermi National Accelerator Laboratory, Batavia, IL 60510 
^Department of Astronomy & Astrophysics, 
The University of Chicago, Chicago, IL 60637 and 
"^Kavli Institute for Cosmological Physics, Chicago, IL 60637 
(Dated: June 14, 2011) 

The Planck Collaboration measured the Sunyaev-Zel'dovich (SZ) decrement of optically selected 
clusters from the Sloan Digital Sky Survey, finding that it falls significantly below expectations 
based on existing mass calibration of the maxBCG galaxy clusters. Resolving this tension requires 
either the data to go up, or the theoretical expectations to come down. Here, we use data from 
the Wilkinson Microwave Anisotropy Probe (WMAP) to perform an independent estimate of the 
SZ decrement of maxBCG clusters. The recovered signal is consistent with that obtained using 
Planck, though with larger error bars due to WMAP's larger beam size and smaller frequency 
range. Nevertheless, this detection serves as an independent confirmation of the magnitude of the 
effect, and demonstrates that the observed discrepancy must be theoretical in origin. 



Introduction. Galaxy clusters are the largest grav- 
itationally bound objects in the Universe and po- 
tentially powerful probes of dark matter and dark 
energy. Understanding clusters well enough to 
use them as probes requires a multi-wavelength 
approach, as observations from the radio to the 
gamma ray end of the spectrum reveal differ- 
ent features of clusters. Observations of the cos- 
mic microwave background (CMB) can detect the 
Sunyaev-Zel'dovich (SZ) effect 01 HIS], wherein 
hot electrons distort the spectrum of passing pho- 
tons. The distortion has a characteristic spectral 
shape and the morphology of the signal traces the 
integrated pressure of the cluster. 

The SZ effect was first detected with high reso- 
lution and sensitivity measurements of individual 
clusters, and these have become so powerful that 
hundreds of clusters are expected to be detected by 
telescopes such as the South Pole Telescope (SPT) 
[5] and Atacama Cosmology Telescope (ACT) [7]. 
However, with large galaxy surveys such as the 
Sloan Digital Sky Survey (SDSS), which detect 
tens of thousands of clusters optically, even rela- 
tively low resolution CMB surveys such as WMAP 
[S] can detect the SZ signal by stacking many clus- 
ters with the same optical properties. 

Recently, the Planck Collaboration [I] per- 
formed such a measurement utilizing the SDSS 
maxBCG cluster catalog [TU]. This catalog has 
been used to place tight cosmological constraints 
on the amplitude of matter fluctuations [15] that 
are consistent with those derived from X-ray se- 
lected cluster catalogs [T7HIS]- Moreover, the cor- 
responding cluster mass estimates from Johnston 
et al. [12 — modified as per Rozo et al. [TS] — are 



consistent with the X-ray luminosity as estimated 
in [14j , and velocity dispersion measurements as es- 
timated by [T3|. Surprisingly, the Planck Collabo- 
ration found that the SZ-decrement from these 
clusters fell significantly below their expectations. 
If real, this discrepancy signals that either some 
aspect of the physics of the intra-cluster medium 
is not properly understood, that the selection of 
optically selected clusters is not adequately under- 
stood despite the successes it has enjoyed thus far, 
or that there is a hidden systematic in the data. 
We note too that this discrepancy may not the 
same as that discussed in other works that com- 
pare the predicted and observed SZ signal of X-ray 
selected clusters (e.g. |S1 I^TH^ . though see also 
[Ml 123), as [IH] exphcitly checked that their SZ 
estimates are in fact consistent with expectations 
from X-ray. 

Here we attempt to confirm the Planck mea- 
surement by searching for the SZ signal in WMAP 
caused by the SDSS maxBCG galaxy clusters. We 
analyze the data using both parametric and non- 
parametric techniques. In both cases, we detect 
a signal that increases with cluster richness, with 
an amplitude consistent with that observed by 
Planck (Ij. Throughout this work, we adopt a 
fiducial flat ACDM cosmology with 0,^ = 0.3 and 
h — 0.7. Cluster masses are defined within a ra- 
dius -R500 such that the overdensity is 500 times 
that of the critical density of the universe. 

Data. From the WMAP single-year, foreground- 
reduced, W-band, high-resolution {Nside — 1024) 
temperature fluctuation maps, we take the noise- 
weighted mean to construct 7-year coadded maps. 
For each year, we subtract off an overall map 
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mean found from the noise-weighted average of 
10^ random, unique, unmasked pixels containing 
no MaxBCG clusters and removing outliers. The 
resulting map contains the temperature and noise 
in 12.6 million (3.4')^ pixels. 

We use the WMAP maps to estimate the SZ sig- 
nal of maxBCG galaxy clusters [10]. The catalog 
is constructed by searching for red-sequence galax- 
ies drawn from the SDSS Data Release 5, and in- 
cludes 13,823 galaxy clusters over 7500 deg^ in a 
redshift slice z € [0.1,0.3]. Each cluster is assigned 
a richness N2Q0, the number of galaxies within 2a 
of the red-sequence and within a specified radial 
aperture [seelTUl for details]. Cluster redshifts are 
photometrically estimated and are accurate at the 
Az < 0.01 level. The input for our analysis is the 
position angle {RA, DEC), redshift (z), and rich- 
ness {N200) of each cluster. 

Analysis. The SZ signal at angular position 6 
relative to a cluster center is an integral of the 
pressure P along the line of sight: 
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where g is the characteristic spectral shape of the 
SZ distortion equal to —2 at very low frequencies, 
ctt is the Thomson cross section, and — with the 
cluster centered at the origin — the 3D position x 
has components equal to {dA{z)9,l), where dA{z) 
is the angular diameter distance out to the redshift 
of the cluster. 

We employ two complementary techniques to ex- 
tract the weak SZ signal from the WMAP data: a 
parametric and a non-parametric approach. In the 
parametric ( "matched filter" ) approach, we take as 
input the data di = AT^ in pixels around each clus- 
ter, the covariance matrix Cij connecting the pixels 
(which includes both instrumental noise and cor- 
relations due to the primordial CMB), and a mor- 
phological template ti for the signal in each pixel. 
We then minimize with 
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where A is the dimensionless amplitude of the sig- 
nal. The sum is over all A'pixcis within an angular 
radius Omax = 5 x ^500 of the cluster center, where 
^500 — Rboo/dAiz). The estimate of -R500 requires 
a mass-richness relation, which we take to be (HH) 
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^500 


^00 X 104 


m) 


Count 


(IQI^Mq) 


(arcmin) 


(arcmin^) 


8-30 (14) 


12408 


8.5 


3.2 


0.33 ± .23 


31-50 (37) 


690 


24 


4.5 


2.4 ± 1.3 


51-80 (60) 


160 


40 


5.3 


4.5 ±2.8 


81-200 (104) 


33 


71 


6.4 


17 ±5.9 



TABLE I. Properties of the four bins in iV2oo. The 
number of clusters is taken after the WMAP mask is 
applied, and the mean number in each bin is used to 
calculate M500 and 61500 for a redshift z — 0.2. The SZ 
measurement is given in the final column. The errors 
are estimated from the filter noise. 



with Bm\n = 0.95 ±0.12 and aA/|Ar = 1.06 ±0.11. 
Table I lists the richness bins in 7V2001 the assumed 
masses, and the angular extent of the clusters at 
z = 0.2. 

To obtain the elements of the covariance ma- 
trix, we first compute the beam-convolvec{^ CMB 
correlation function ^{6) using the C/ correspond- 
ing to WMAP 7 parameters [6]. For each pair 
of pixels z, j we determine the angular separation 
= [ef + 6*1 - 20,0 j cos(0, - 0j)]^/^ between the 
centers of the pixels, and set the element of the 
covariance matrix Cij — £,{9ij)- The full covari- 
ance matrix is set by adding the diagonal noise 
contribution NfSij. For the template, we choose 



the pressure profile in Arnaud et al 
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where p{x) is the dimensionless profile 

8.403 
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with C500 = 1.177, a = 1.051,7 = 0.3081, and 
6 ~ 4.931. The predicted signal in a pixel i is 
ATsz as computed from Eq. ([T]), convolved with 
the WMAP beam B{0-0,): 
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Minimization of the for a given cluster a leads 



r, ( iV^nn \ , . 1 We take a Gaussian beam of variance (0.1°)'^ for the W 
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frequency band. 



3 



to an estimate of the amplitude: 

ij 

with "fiher noise" (the error due solely to instru- 
ment noise+CMB) defined as 



An amplitude estimate Aa — 1 corresponds to data 
that is exactly consistent with the assumed tem- 
plate. The mean amplitude (A) for a richness bin 
is defined as the inverse- variance weighted average 
of Aa over all clusters in the bin. Due to intrinsic 
scatter in the mass-richness relation, the total error 
on A will be larger than the filter noise. Therefore, 
we inverse-variance weight the amplitudes by the 
total error, defined as the sum in quadrature of the 
filter error and the intrinsic scatter. To estimate 
the intrinsic scatter, we compute the scatter in a 
given richness bin. 
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and then subtract cr^^^ juter- Since the weighting 
depends on the total error, we iterate until the 
input total error is equal to the output scatter. 

The template and the measured value for the 
amplitude A can be transformed into an estimate 
of any related quantity. One relevant quantity is 
5^500) defined as 
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where the integral is restricted to r < -Rsoo- 
the template in Eq. ([2]), we find 



Using 



F500 = 8.5 X IQ-^A 



^ 3 X 1014Mq 



1.79 



arcmm 



We also carry out a non-parametric fit by aver- 
aging the temperatures in pixels binned into an- 
nuli around the clusters. As in the matched fil- 
ter approach, the covariance matrix for the annuli- 
averaged temperatures contains contributions from 
instrument noise and from the primordial CMB. 
The first is straightforward to propagate using the 
combined weights from the single year WMAP 
maps. To compute the second, we first note that 
the estimated temperature in a given annulus a 
distance 9 from a given point, taken to be the ori- 
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FIG. 1. W band SZ signal in 4 iVaoo bins. The bine 
stars show the expectation for the SZ decrement of 
the maxBCG clnsters based on the Rozo et al. (2009) 
mass-richness relation, while the red circles show the 
signal measnred nsing WMAP data. The error bars are 
set by the filter noise. The solid line shows the best-fit 
result from Planck. The inset shows that randomiz- 
ing the positions of the clusters within the survey area 
leads to a signal consistent with zero. 



gin, isTg = SiLi T{9, (^i), where the sum is over 
all N pixels in the annulus, each identified by its 
polar position (pi. The covariance matrix between 
measurements in two annuli at 9 and 9' is 

N N' 



Here N' is generally not equal to N since larger 
annuli have more pixels (of fixed size) in them. 
The average of the product of the temperatures is 
£,{9ii'), the correlation of the CMB between two 
pixels separated by angular distance 9iii . 

We use 11 annuli with angular distance from the 
cluster equal to (0.5, 0.9, 1.3, . . . , 4.5) x 6*500. The 
CMB noise is the same order of magnitude as the 
instrumental noise but is highly correlated between 
rings, so a simple ^-by-eye of the profile is not suf- 
ficient to interpret the results. We ran 100 mocks 
to test this approach and found that the extracted 
amplitude A is biased high by 1.35 ± 0.25 for the 
51 — 80 richness bin and 1.67±0.62 for the top bin. 
We correct for these biases below. 
Results. Fig. [l] shows our estimates for 1500 for 
clusters in the four richness bins. The corre- 
sponding amplitudes are A = 0.37 ± 0.26, 0.42 ± 
0.23,0.32 ± 0.20,0.43 ± 0.15. The error bars are 
(J (A). filter propagated from A to Y^qq. We also esti- 
mate the total scatter with (T(^A),totah which is sen- 
sitive to outlier clipping. In Table [III we values 



for 



^ (A), filter 



and (T(^ji-^ fatal without clipping, and 



compute the contribution from intrinsic scatter. 



4 
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^{A) ^intrinsic 


14 


0.60 


0.23 


0.55 


37 


1.8 


1.3 


1.2 


60 


4.3 


2.8 


3.2 


104 


11 


5.9 


9.3 



TABLE II. Mean Y^oq errors in units of 10"'* arcmin^. 
The filter noise includes only uncertainties from the 
instrument and the CMB. The quadrature difference 
between the total dispersion and the filter noise gives 
an estimate for the error in the mean due to intrinsic 
scatter. 




Our measurements are in excellent agreement 
with those of the Planck Collaboration |Tj, and 
are clearly lower than the theoretical predictions. 
Given this independent confirmation of the ten- 
sion first observed in [1], it is clear that either the 
intra-cluster properties of galaxy clusters are not 
adequately understood, or that the optical cluster 
selection suffers from a systematic that has gone 
undetected. We do not currently know the solu- 
tion to this problem, and indeed, the discrepancy 
is likely to remain an active field of research in 
the immediate future. We note that a similar dis- 
crepancy appears to arise for the SZ signal about 
Luminous Red Galaxies (LRGs) [2U] . 

We perform two systematics cross-checks of our 
analysis pipeline: first, we randomize the location 
of every cluster in the catalog and repeat our mea- 
surement. As expected, we find no signal. In addi- 
tion, we generate Monte Carlo realizations of the 
experiment to test whether the recovered ampli- 
tude is biased. We randomly sample the maxBCG 
catalog, and assign to every cluster an SZ decre- 
ment (including scatter, for which we vary crinY 
from 0.6 to 1.2). The cluster signal of Eq. 
is convolved with the beam and added to ran- 
domly generated (beam-convolved) CMB maps. 
We then fit for the amplitude of the signal using 
our pipeline, finding that the recovered signal is 
biased low by 2% (3%) in the 51 - 80 bin and 5% 
(6%) in the 81 - 200 bin for ai^Y = 0.6 (1.2). 

Fig. [2] shows the non-parametric measurement 
for the signal as a function of distance from the 
cluster center in the two largest richness bins. A 
simple assuming no signal reveals strong detec- 
tions in both bins [x^ = 22.3, 24.2 for 10 degrees 
of freedom each). The smooth curves show the 
predicted signal from our fiducial template tem- 
plate {A = 1). The best fit for the two bins yields 
A = 0.67 ± 0.22 for the Ngais = 51-80 bin and 
A = 0.67 ± 0.20 for the highest richness bin. Cor- 



FIG. 2. W band SZ signal as a function of distance 
from the cluster center for clusters in the two largest 
richness bins listed in Table 1. Error bars include in- 
strumental noise only. Solid curves are the predictions 
from the template assumed in Eq. ([2|. The deviation 
on small angular scales translates into a smaller am- 
plitude than assumed, in agreement with the Planck 
results. 



reefing for the bias in these bins leads to estimates 
of A = 0.50 and 0.40 respectively, consistent within 
errors with the parametric determination. The lat- 
ter is the superior way to determine the amplitude 
as it weights all pixels optimally, albeit at the cost 
of being sensitive to the assumed template. 

SD and JH are supported by the US Department 
of Energy, including grant DE-FG02-95ER40896; 
and SD by National Science Foundation Grant 
AST- 0908072. ER is funded by NASA through the 
Einstein Fellowship Program, grant PF9- 00068. 



[1] N. Aghanim et al. [Planck Collaboration], 

arXiv:1101.2027. 
[2] M. Arnaud et al. 2009 arXiv:0910.1234. 
[3] Birkinshaw, M. 1999, Phys. Reports, 310, 97 
[4] Carlstrom, J. E., Holder, G. P., & Reese, E. D. 

2002, Ann.Rev.Ast.Astrophy., 40, 643 
[5] N. Jarosik et al, Astrophys. J. Suppl. 192, 14 

(2011). 

[6] E. Komatsu et al. [WMAP Collaboration], Astro- 
phys. J. Suppl. 192, 18 (2011). 

[7] F. Menanteau et al, Astrophys. J. 723, 1523 
(2010). 

[8] Sunyaev, R. A. & Zeldovich, Y. B. 1972, Com- 
ments on Astrophysics and Space Physics, 4, 173 

[9] K. Vanderlinde et a/., Astrophys. J. 722, 1180 
(2010) 

[10] Koester, B. P. et al. 2007, Astrophys. J., 660, 239 
[11] Sheldon, E. S. et al. 2009, Astrophys. J., 703, 2217 
[12] Johnston, D. E. et al. 2007, arXiv:0709.1159 
[13] Becker, M. R. et al. 2007, Astrophys. J., 669, 905 



[14] RykofF, E. S. et al. 2008, Astrophys. J., 675, 1106 
[15] Rozo, E. et al. 2010, Astrophys. J., 708, 645. 
[16] Rozo, E. et al. 2009, Astrophys. J., 699, 768. 
[17] Mantz, A. et al. 2010, MNRAS, 406, 1759. 
[18] Vikhlinin, A. et al. 2009, Astrophys. J., 692, 1060. 
[19] Henry, J. P. et al. 2009, Astrophys. J., 691, 1307. 
[20] Hand, N. et al. 2011, arXiv:1101.1951. 
[21] Lieu, R. et al. 2006, Astrophys. J., 648, 176. 



5 



[22] Biolby, R.M. and Shanks, T. 2007, MNRAS, 382, 
1196. 

[23] Diego, ,J.M. and Partridge, B. 2010, MNRAS, 402, 
1179. 

[24] Afshordi, N. et al. 2007, MNRAS, 378, 293. 
[25] Mehn, J.B. et al. 2011, A&A, 525, 139. 
[26] Planck Collaboration, arXiv:1101.2043. 



